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We consider the necessary conditions of superconducting pairing at repulsive interaction between 
particles composing a pair with large total momentum: (1) the existence of, at least, one negative 
eigenvalue of the repulsion potential and (2) mirror nesting of the Fermi contour. Under these 
conditions, we represent the solution of the self-consistency equation continuously depending on the 
momentum of the relative motion of the pair. The corresponding superconducting order parameter 
changes its sign on a line crossing the Fermi contour inside the domain of definition of the relative 
motion pair momentum. We argue that repulsive-induced superconducting pairing with large total 
pair momentum may be just the case relating to high-temperature superconducting cuprates. 

PACS numbers: 78.47.+p, 78.66.-w 

Keywords: superconductivity, repulsion interaction, nesting, large pair momentum, cuprates 



I. INTRODUCTION 

A conventional point of view relating to the mecha- 
nism of high-temperature superconductivity of cuprates 
is absent up to this point. Nevertheless, there is an exper- 
imental evidence that there is the singlet pairing of cur- 
rent carriers below the superconducting (SC) transition 
temperature in these compounds although the question 
is what interaction dominates pairingi In this connec- 
tion, the models of the superconductivity both with dom- 
inating attraction and repulsion are considered^ Another 
open question is what contribution into the pair binding 
energy (due to potential energy as it is in the conven- 
tional theory by Bardeen, Cooper and Schrieffer (BCS) 3 
or, maybe, due to kinetic energy^) may be considered as 
a driven force of the SC transition in the cupratest^ 

In contrast to the conventional paradigm that the SC 
transition is a consequence of an instability of the Fermi 
liquid with respect to arbitrarily weak attraction between 
particles the question relating to the character of the 
ground state of normal (N) phase of cuprate compounds 
remains under discussion^ The manifestation of the 
pseudogap in underdoped cuprates is a reason to suppose 
that the N phase may be something different from the 
normal Fermi liquid. Assuming that cuprates are strong 
correlated electron systems one can, in principle, connect 
rather wide spectrum of physical properties of these com- 
pounds which include, besides typical of cuprates super- 
conductivity and antiferromagnetism (AF), some other 
ordered states arising due to a competition of SC and 
AF orders^ In this connection one can raise the ques- 
tion of the applicability of the mean-field theory (like the 
BCS theory which describes conventional superconduc- 
tors successfully) to cuprate superconductors A 

Pairing with zero total pair momentum does not im- 
pose any kinematical constraint on the momenta of the 
relative motion coinciding, in this case, with the mo- 
menta of the particles composing the pair. In the BCS 
model with an attractive interaction between particles, 
however, there is a "dynamical" constraint connected 



with the fact that the domain of the effective attraction 
is a narrow layer enveloping the Fermi surface. The vol- 
ume of this layer is proportional to the statistical weight 
of the pair that is a number of one-particle states which 
contribute to the state of the relative motion of the pair. 9 

On the contrary, at non-zero total pair momentum 
there is an essential kinematical constraint resulting in a 
finite domain of the momentum space in which momenta 
of the relative motion of the pair should be defined^ A 
supplementary dynamical constraint arising from an at- 
tractive interaction, decreases such a domain abruptly re- 
sulting in a considerable decrease in the statistical weight 
of the pair. It should be emphasized that, in the case of 
Coulomb repulsion, any dynamical constraint is absent 
and the statistical weight of the pair is fully determined 
by the volume of kinematically allowed domain of the 
momentum space. 

The Fermi surface of quasi-two-dimensional (2D) elec- 
tron system of copper-oxigen planes of a cuprate com- 
pound is degenerated into a line which we call the 
Fermi contour (FC). In cuprates, the FC is situated in 
an extended vicinity of a saddle point of the electron 
dispersion. 11 As a consequence, one may expect a FC 
with strong nesting feature and, correspondingly, a con- 
siderable difference in the transversal (with respect to the 
FC) and the longitudinal Fermi velocities^ In addition, 
one can assume that the corresponding effective masses 
should be of opposite signs and differ in absolute values 
considerably. Thus, one can expect that, for some def- 
inite (antinodal) directions coinciding with the sides of 
2D Brillouin zone, there exists rather large area of 2D 
momentum space defining the statistical weight of the 
pair with certain non-zero total momentum^ One can 
create a current-less state corresponding to such a pair 
as a definite linear combination of the states relating to 
crystal equivalent momenta^ 

The SC instability of 2D electron system arises due 
to the fact that the boundary separating occupied and 
unoccupied one-particle states is a line. In the case of 
zero total pair momentum, this line is the full FC whereas 
at non-zero total pair momentum such a boundary is, 
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generally speaking, a set of points. Therefore, in such 
a case, there is no logarithmic singularity in the self- 
consistency equation and, consequently, the SC pairing 
is impossible at arbitrarily weak interaction strength. 

However, the extended saddle point may result in a rise 
of such a topology of isolincs in a vicinity of the FC that 
the boundary between occupied and unoccupied states 
of the pair relative motion turns out to be a set of lines 
of finite length. These lines form a "pair" Fermi contour 
(PFC) for the relative motion of the pair with definite 
total momentum. Such a feature (mirror nesting) of the 
FC— is held in a finite doping interval and may warrant 
a logarithmic singularity in the self-consistency equation 
resulting in the SC instability at arbitrarily interaction 
strength. 

The channel of the SC pairing with large total pair mo- 
mentum results quite naturally from a consideration of 
the competition of Cooper (at zero pair momentum) SC 
pairing and AF ordering in the framework both a band 
modeii^ and the models which concern strong electron 
correlations. 16 Also, such a channel arises with neces- 
sity in rather general phemenology based on symmetry 
considerations— 

Recent experimental results relating to neutron scat- 
tering in cuprates (the so-called 41 meV peak) may be 
considered as an indirect evidence in favor of the SC 
pairing with large total momentum— Another two fun- 
damental experimental results, the so-called "peak-dip- 
hump structure" of angle-resolved photoemission spec- 
trum (ARPES) 1 ^^ and relatively small superfluid phase 
stiffness typical of cuprates^ can be explaned qualita- 
tively in the framework of the conception of mirror nest- 
ing and repulsion-induced SC pairing with large total pair 
momentum as well. 

The goal of this paper is an investigation of the neces- 
sary conditions of the repulsion-induced SC pairing with 
large total pair momentum and a solution of the cor- 
responding self-consistency equation arising within the 
mean-field scheme. 



II. MIRROR NESTING 

The kincmatical constraint due to the presence of the 
FC results in the fact that both momenta fc ± of like- 
charged particles composing a pair with total momen- 
tum K = k + + k have to be situated either inside or 
outside the FC. Thus they belong to dependent on K 
domain of the momentum space which may be consid- 
ered as a statistical weight of the pair. This domain R K 
is symmetrical with respect to the inversion transforma- 
tion k <-» — k of the momentum of the relative motion of 
the pair, k = (k , — k_)/2, and, generally speaking, con- 
sists of the two parts sfc"^ and inside and outside 
the FC, respectively. An excitation energy of the pair 
with respect to the chemical potential \x that is kinetic 



energy of the particles composing the pair, 

2£ Kk = e(K/2 + k) + e(K/2-k)-2», (1) 

equals zero on the boundary separating occupied (S^- ^ ) 

and unoccupied (Sj^ ) parts of the domain S^-. Here- 
after, for the sake of simplicity, we restrict ourselves 
to a consideration of hole-doped cuprate superconduc- 
tors within a simple single-band scheme; therefore, s{k±) 
makes sense of a hole dispersion. 

As it is mentioned above, at K = 0, the momenta of 
the relative motion ±fc coincide with the momenta of the 
holes composing the pair and therefore the boundary sep- 
arating the occupied and unoccupied parts of the domain 
S K is the full FC whereas at any K ^ the boundary 
between S^- ^ and S»t* , generally speaking, turns out to 
be a set of isolated points. However, due to a special 
hole dispersion, such a boundary, at definite total pair 
momentum, may coincide with finite-length pieces of the 
FC as it is shown schematically in Fig.l. In such a case 
the boundary between S^- ^ and plays role of a pe- 
culiar FC with respect to the relative motion of the pair 
and may be called a pair Fermi contour (PFC). 

Using the symmetry property, k — fc, one can write 
the condition defining the PFC in the form 

e{K/2 + k) - [i = e{K/2 - k) - fi, (2) 

under the additional condition that both k and —k be- 
long to the domain 2^. It is obvious that, to obey the 
condition J2Jl , the vector K should be directed along one 
of the symmetrical (antinodal or nodal) direction of 2D 
Brillouin zone. Hence the boundary separating occupied 
and vacant states has to obey a mirror symmetry with 
respect to the pair momentum K direction and so Eq. J2J) 
may be called a mirror nesting condition. It should be 
also noted that, under mirror nesting condition, the in- 
version transformation k «-> k superposes an occupied 
state with an occupied one and a vacant state with a 
vacant one, respectively. 

In the framework of a single-band scheme relating to 
cuprates, one can give some rather obvious examples of 
the mirror nesting condition resulting in a rise of the 
PFC— i2i In a case when the FC of a hole-doped cuprate 
compound is simply connected and appears as a square 
with rounded corners^i one can expect that a change 
of FC curvature sign, which can be easily realized be- 
cause of the existence of long, almost rectilinear, pieces 
of the FC, may provide a realization of the mirror nesting 
condition for certain chosen pair momenta. 21 One more 
single-band example represents an underdoped cuprate 
compound considering as doped Mott insulator in which 
the FC encloses some symmetrically disposed hole pock- 
ets around the points belonging to symmetrical directions 
in 2D Brillouin zone— In such a case, the mirror nest- 
ing condition is perfectly satisfied when a half of the total 
pair momentum corresponds to the center of a pocket and 
the whole of the line enclosing the pocket turns out to 
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FIG. 1: Domain of definition S K — "S K '+5^ of pair relative 
motion momentum (schematically; it is shown only a half of 
the domain): the subdomain a K ' (H^"') corresponds to the 
shaded occupied (unoccupied) part of the domain S. These 
subdomains are separated by the PFC (solid line) being a part 
of the FC (fine line). The nodal line k — k (dashed line) 
separates the parts of the momentum space corresponding to 
different signs of the order parameter. 



be the PFC. The stripe structure arising in underdoped 
cuprates gives one more possibility resulting in a rise of 

p FC 20 



III. THE SELF-CONSISTENCY EQUATION 

To consider the problem of pairing with non-zero total 
pair momentum in the mean-field framework let us write 
the Hamiltonian corresponding to pair relative motion in 
the form 



H 



k 



i a k'-i a k' + ] 



(3) 



where e k± = e(k±), U(k — k') is the Fourier transform 
of the interaction energy, S is a normalizing area, a,/, 
(a k±a .) creates (annihilates) a hole with the momentum 
k± and the spin quantum number a — t, I- The sum- 
mation in Eq.© is taken over all range of the allowed 
values of the momentum of the relative motion and thus 
is restricted by the domain S/f . Note that the summa- 
tion in the Hamiltonian Eq. 10 should be taken over only 
two (instead of three in a general case) variables, k and 



k', just as in the case of BCS Hamiltonian, taking into 
account the only interaction between particles composing 
the pairs with total pair momentum K. 

As usual, to diagonalize the Hamiltonian Eq.© ap- 
proximately one can introduce creation and annihila- 
tion operators of new quasiparticles using the well-known 
Bogoliubov-Valatin transformation: 



a fe+T 



l Kk b k, + l 



Kk u k,-1 



Kk°k,+V 



V, 



(4) 



The Hamiltonian, up to the terms of the order of b 2 , can 
be written as 



H 



(0) 



H 



(i) 



The ground state energy has the form 



E 



AO 



% ^ £,Kk v Kk 
k 



k 



A 



Kk Ul Kk u Kki 



(5) 



(6) 



where, related to the value of the chemical potential, an 
energy of the relative motion of the pair with the total 
momentum K is defined by the Eq.© and the order 
parameter is defined as 

-^U(k-k')u Kk ,v Kk ,. (7) 



A' A- 



k' 



Diagonal and nondiagonal, with respect to quasiparticlc 
operators, parts of the Hamiltonian can be written as 



Mo) _ 

K ~ 



k-J3=±l 



(8) 



and 



o-(l) _ 



E 

k 



[2* 



Kk a Kk u Kk 



( v Kk 



l Kk 



)A 



AfeJ 



x ( b k,+i b k,-i+h,-ib K +i) 



(9) 



respectively. Here the energies corresponding to two 
branches (/3 = ±1 ) of the one-particle excitation spec- 
trum are equal to each other, 



A 2 



(10) 



A choice of the amplitudes, u Kk and v Kk , in 
Bogoliubov-Valatin transformation Eq. Q is determined 
in the zero-temperature limit by the conditions that (i) 
all of the states inside the subdomain , in which the 
kinetic energy of the relative motion of the pair is nega- 
tive, 2£ff fc < 0, must be occupied and (ii) the nondiagonal 
part, Eq.©, of the Hamiltonian vanishes. In addition, 
the condition u 2 Kk + v 2 Kk = 1 preserving Fermi's com- 
mutation relations for quasiparticle operators must be 
fulfilled. These conditions yield 



1 



u Kk 



1 - 



A A- 



V %Kk + ^Kk , 

1 A, 
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(11) 
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It should be noted that the Bogoliubov-Valatin ampli- 
tudes 1(11(1 are determined not by the one-particle energy 
e(k) — fj, as it were at K — but by the relative motion 
energy Q which transforms into s(k) — fi just in the case 
K = 0. 

In the zero-temperature limit, one can obtain the self- 
consistency equation determining the SC order parame- 
ter, 



x Kk 



-y 

2S ^ 



U(k-k')A 
7F 



Kk' 



Kk' 

"A 2 " 

^Kk' 



(12) 



where k and k' are momenta of the relative motion of 
pairs with one and the same total momentum K, U(k) a 
matrix element of the effective interaction potential in the 
Hamiltonian Eq.(PJ depending on the momentum trans- 
fer k — k — k' due to a scattering from an initial (k) to 
a finite (fc') state of the relative motion of the pair being 
also a difference, k + — k' + , between the corresponding 
momenta of the particles composing pairs with the total 
momentum K . Summation over momenta of the rela- 
tive motion in the equation l|12|l is performed within the 
whole of their domain of definition S^. Thus, the SC 
order parameter is defined inside the domain S^- as well. 

To solve the equation (|12|) it is convenient to reduce 
it to the corresponding integral equation considering the 
quantities A Kk — A(fc) and ^Kk = £,(k) as continuous 
functions of the relative motion momentum. Hereafter, 
for a simplification of all of the notations, we omit the 
label K . Using the notations 



V (k) = V£ 2 (*) + A2(fc), /(fc) = l/47r7j(fc), 
one can rewrite ()12|) in the form 

A(fc) = -^f u (\ k - fe'|)/(fc')A(fc')rf 2 fc', 



(13) 



(14) 



where U(\k — k'\) may be called a kernel of the non-linear 
integral equation ljT4*|> . 



IV. NECESSARY CONDITIONS OF 
REPULSION-INDUCED SC PAIRING 

In this section, we discuss the necessary conditions of 
the repulsion-induced SC pairing with large total pair 
momentum and put forward a procedure of reducing of 
the mean-field SC gap equation 1(141) to an approximate 
non-linear singular integral equation with a degenerate 
kernel^ 

The matrix element U(k — fc') is connected with the 
effective interaction potential U{r) in the real space by 
the Fourier transform 



U(\k-k'\ 



U(r) exp [i(k 



k')r]d 2 i 



(15) 



One can deduce a simple criterion of the absence of non- 
trivial solutions of the equation i|14|) . Let us multiply 



this equation and the function /(fc)A(fc) together and 
integrate the obtained relation over fc inside the domain 
H. Expressing U(k — fc') in accordance with 115fl and 
changing the order of integration over fc, fc' and r one 
can obtain 



f A 2 (k)f(k)d 2 k = -i- J U{r)L{r)d 2 r, 



where 



L(r) 



J A(k)f(k)exp{ikr)d 2 k 



(16) 



(17) 



is a nonnegative function of r. The integration over r 
in (|16(l is performed over the whole of 2D real space. As 
far as /(fc) > 0, the left side of (|16|l is nonnegative as 
well. If the potential U(r) were positive at any r, the 
right-hand side of the Eq. (|16fl . on the contrary, turns 
out to be negative and the equality 1(16(1 may be satisfied 
only in the trivial case when A(fc) = in the whole 
of the domain 2. Thus, under the condition that the 
interaction is purely repulsive in the real space that is 
U(r) > at any r, the self-consistency equation leads 
only to the trivial solution. 

Therefore, a rise of repulsion-induced SC order is pos- 
sible only in the systems where the interaction potential 
is a function of r with alternating signs. As an example of 
such a potential with alternating signs, one can consider 
screened repulsive Coulomb interaction in a degenerated 
Fermi system which exhibits the Friedel oscillations. 24 
Further, we restrict ourselves to a consideration of such 
a repulsive potential U(r). 

One can see that the maximal (positive) value of the 
function r[/(r), arising when the integral in the right- 
hand side of 1)16(1 is calculated using the polar coordi- 
nates, corresponds to r = 0. This function with alter- 
nating signs is multiplied by a nonnegative function L{r) 
defined by Eg. 1(17(1 . The governing contribution into the 
product A(fc)/(fc) determining the function L(r) is due 
to a small vicinity of the PFC because of the equality 
£(fc) = which takes place just on the PFC. 

As it will be shown below, the gap function A(fc) re- 
sulting from the repulsion interaction changes sign on a 
line crossing the PFC. Therefore, the absolute value of 
the integral in ((17(1 represents a function of r appear- 
ing as a series of peaks with decreasing magnitude and 
a character separation between the neighboring peaks of 
about ir/kp where k P is of about the length of the PFC. 
The first (and the largest) peak results in the fact that 
the product of the factor A(fc)/(fc) with alternating signs 
and exponential function in the integral 1(17(1 turns out to 
be a function with the values mainly of constant sign. 

The value L(0) is, generally speaking, a small quan- 
tity because of a partial (or even full) compensation of 
the contributions of the parts of the domain S in which 
A(fc) has opposite signs. For example, in the case of the 
full compensation, L{r) ~ r 4 , therefore the main pos- 
itive extremum of the function rU (r) turns out to be 
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FIG. 2: Behavior of the functions rll(r) and L(r) contribut- 
ing to the right-hand side integral in Eq. ll6ll . schematically. 
One can see L(r) as the function l|17fl averaging over the angle 
variable. 



suppressed considerably. In the opposite case of large r, 
the function L(r) is also small due to rapidly oscillating 
factor exp(ikr). As far as k P is pronouncedly less than 
k F , the position of the first maximum of the function 
Q17[l. r 1 oj it /k p > n/kp, corresponds quite naturally 
to the region of the Friedel oscillations of the repulsive 
Coulomb potential. There is a significant region of the 
r - space in which the integrand in the right hand side of 
the Eq. is negative whereas the contributions of the 
regions with positive integrand are comparatively small 
as it is illusyrated schematically in Fig. 2. Therefore the 
integral in may turn out to be negative. 

In such a case, the equality (|16|l is provided due to the 
fact that the function /(fe) is singular on the PFC and 
may be arbitrarily large at A(fc) — » 0. The magnitude of 
the oscillations of the function L(r) depends on A(fc) be- 
cause an effective width of the vicinity of the PFC leading 
to the main contribution into the integral (|16|l is deter- 
mined by just the function A(fc). Thus, the condition 
may be fulfilled at sufficiently small A(fe) ^ 0. 

Therefore the presence of the PFC which results in the 
singularity of the Eq. I|14fl is the first necessary condition 
of the repulsion-induced SC pairing. 

The equation l|14|l is a non-linear integral Hammerstein 
equation with a symmetric kernel U(\k — k'\). All of the 
eigenvalues of such a nondegenerate kernel are real and 
belong to a discrete spectrum \ n (where n = 1,2,...) 
with the condensation point |A n | — > oo at n — > oo. If a 
kernel U(\k — k'\) is the Fourier transform (|15|) of every- 
where positive effective potential U (r) (a positive defined 
kernel) all its eigenvalues are positive as well. In such a 
case, as it follows from Ea. (|16fl . the equation l|14|) has the 
trivial solution only. 



Therefore, the presence of at least one negative eigen- 
value in the spectrum of the kernel U(\k — k'\) is the sec- 
ond necessary condition of the existence of a non-trivial 
solution of the equation (|14|) . 

It should be noted that, as it follows from Eci. (|14fl . 
a non-trivial solution A(fe) ^ of the self-consistency 
equation at U{k — k') > 0, if such a solution exists, must 
be a function with alternating signs inside its domain of 
definition S. 

In the case of the non-degenerate kernel U(\k — k'\), 
the existence of such a non-trivial solution of the self- 
consistency equation at the repulsive interaction can be 
demonstrated if one represents this kernel as an expan- 
sion over its eigenfunctions <p n (k) which are the solutions 
of the linear integral equation 

¥>„(*) = K l U (\ k - k'\)Vn{k')d 2 k'. (18) 

Let us write this expansion in the form 25 

u{k _ k>) = £ rt(fc)¥>n(fc') (19) 

n ™ 

and substitute it into the self-consistency equation H14[) . 
Integrating the right-hand side of this equation term-by- 
term, one can see that a non-trivial solution, if it exists 
at all, may be represented in the form of an expansion 
over the complete system of the eigenfunctions ip n (k), 

A(fc)=^A„^ n ,(fc), (20) 

n' 

where A n are the expansion coefficients which can be de- 
termined as the solutions of the infinite equation system, 

a ™ = -?tE^'V (21) 

Here, f nn , are the matrix elements of the form 

fnn> = JynWfWtpAk) (22) 

where the function f(k) is determined by Ea. (|13[) . 

Multiplying Eq.J2U and A,jA n together, one can, after 
a summation over n, obtain the relation 

A„A 2 = -i- £ A n / n „,A n , (23) 

n n,n' 

which is analogous to the relation (flfif) . As it is follows 
from Ea. (|23|l . there is only the trivial solution of this 
equation in the case when all of the eigenvalues are pos- 
itive. Therefore, a non-trivial solution may, in principle, 
exist only in the case when at least one of the eigenvalues 
turns out to be negative. In fact, such a condition turns 
out to be not only necessary but the sufficient one. 

To prove this statement, one can take advantage of the 
obvious analogy between the spectral expansion of the or- 
der parameter, Ea. (|20|l . and its expansion into a series of 
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spherical harmonics. As it is well known, the SC pairing 
arises when the negative scattering length corresponds at 
least to one of the values of orbital angular momentum. 26 
For this reason, taking account of such an analogy, we do 
not present a direct proof of such a statement within the 
framework of the spectral expansion. 27 



V. SIMPLE DEGENERATE KERNEL 

The function defining a non-linear operator in Eq. 1141) 
has a special form A/(A) where /(A) may be called a 
nonlinearity factor. In this sense, the equation i|14|) may 
be referred to the class of quasi-linear integral equations. 
In the case of a degenerate kernel, the significant feature 
of such equations is that a form of a non-trivial solution 
(if it exists at all) reproduces the kernel structure. The 
kernel U(\k — k'\) can be approximately reduced to a de- 
generate one if one takes into account that both variables, 
k and k! , are defined in a relatively small domain 2 of the 
momentum space. Therefore, if the function i|15|) can be 
expanded into the Tailor series, one can restrict himself 
to some lowest powers of the argument k = \k — k'\. As 
it follows from the definition (|15fl. the expansion of U{k) 
includes only even powers of k and the Fourier transform 
of the interaction potential can be written as 



U(k) = 2tt 



1 2 , ^ 4 

2 2 T g 4 



where 



U(r) 



- n+1 dr. 



(24) 



(25) 



A degenerate kernel corresponds to a finite number of the 
terms in the expansion l|24|) and may be considered as a 
good approximation of the true nondegenerate kernel in 
the case of a small domain 3. It should be noted that 
the Fourier transform U (n) corresponding to a repulsive 
interaction is, at small decreasing function of its 

argument resulting in u Q > and u 2 > 0. 

To solve the equation (|14|l it is convenient to transform 
it into a dimensionless form. Using two constants U and 
r , being characteristic energy and length scales, respec- 
tively, one can define them by the relations UqVq = u 
and UqTq = u 2 . In the case of screened Coulomb po- 
tential the parameters r and U = e 2 /r make sense of 
a screening length and a characteristic Coulomb energy 
respectively. Thus, the energy U determines the scale 
of the quantities A(fc) and £(fc) whereas a momentum is 
measured by the units of r^ 1 . The equation (|14|) remains 
invariant with respect to such a scaling if one reduces the 
domain of integration 5 to dimensionless variables (the 
components k 1 and k 2 of the momentum of the relative 
motion of the pair). 

If the domain 5 is small enough one may keep only 
two terms of the lowest order in the expansion (|24|l and 
obtain a degenerate kernel in the form 



One can make sure that the obtained degenerate kernel 
has four eigenvalues three of which are positive whereas 
the fourth one is negative. 

To show this one has to write down the linear integral 
equation determining the eigenvalues A and the eigen- 
functions "^(fc) corresponding to the degenerate kernel 
Eq.®, 

tp x (k) = A^ U d (\k - k>\)<p x (k>)d 2 k>. (27) 

It follows immediately from the Eq. I|27|) and the form of 
the kernel l|2t)|l that eigenfunctions must be taken in the 
form 



ip x (k)=a+{ X k)+bk 2 



(28) 



thus reproducing the form of the degenerate kernel l|26|l . 
Here, a and b are scalar coefficients and x is a constant 
vector being subjects to be determined for each of the 
eigenvalues. To obtain the parameters a, b and x of the 
solution (|28|l one has to substitute this solution unto the 
Ea. (l27|) and, after integration, set the coefficients corre- 
sponding to one and the same power of the vector k equal 
to each other. 

It is obvious that the approximation l|2t)|) is sufficient 
if (k + k 1 ) 2 < 2. Therefore, one has to assume that the 
domain 3 is such that the condition 



k 2 < 1/2. 



(29) 



is fulfilled for any momentum of the relative motion. If 
one takes into account that a screening length r is of 
about a few interatomic distances, the restriction l|29|) 
results in the fact that a characteristic size of the domain 
5 should be much less as compared with a characteristic 
Brillouin zone size. 

Due to a symmetry of the domain of integration S with 
respect to inversion transformation k <-» k eigenfunc- 
tions <p^(fc) must be either even or odd. Odd eigenfunc- 
tions are determined by single unknown vector x being 
the solution of the equation 



X = 2n\ {xk')k'd 2 k'. 



(30) 



This equation has two non-trivial solutions directed along 
the symmetry axes of the domain 3. Corresponding 
eigenvalues, 



A,- = 



2ir / kfd 2 k 



i = 1,2 



(31) 



are positive. Here, k i are the components of the vector 
fc along the symmetry axes of the domain H. 

In the case of the even eigenfunctions, one can obtain 
the system of linear homogeneous equations 



2ttA / (l~k' 2 /2)(a + bk iz )d 2 k', 



U d {K) =2n[l-K 2 /2}. 



(26) 



b = ttA / (a + bk' 2 )d 2 k' 



(32) 
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determining two unknown quantities, a and b. The con- 
dition of the non-trivial compatibility of this system re- 
sults in another two eigenvalues of the degenerate kernel 
Eq.JUJl, 



1 



■i -i 



(K - K x ) ± J (K - K,f + (K K 2 - K{) 



Here we use the notation 

K„ I k 2n d 2 k, n= 0,1,2. 



(33) 
(34) 



As it follows from the Cauchy-Schwarz-Bunyakovsky 
inequality for the integrals l|34[l . K K 2 > K\, therefore 
one of the eigenvalues is positive, A + > 0, whereas 
the second one is negative, A_ < 0. Here we take into 
account the fact that, due to the inequality (|29[1 . 2K X < 
K . 

It should be noted that two of the positive eigenvalues 
A 12 ~ "o^ 1 ! another one positive eigenvalue A + ~ u^ 1 , 
and the negative eigenvalue A_ ~ — u Q /u 2 . 

One may neglect the next terms of the expansion l|24() 
only under the condition that, inside the domain H, the 
mean square of the neglected terms is much less in com- 
parison with the absolute value of any eigenvalue of the 
degenerate kernel: 25 




d 2 k I <Pk'{k - k'Y < | AT 



(35) 



In the case of the positive eigenvalues, this inequality is 
satisfied under rather natural condition that the third 
term in the expansion (|24[l is small in comparison with 
the second one. In the case of the negative eigenvalue, the 



inequality (|35|l reduces to the condition that u |u 4 | < cu 2 
where c<Cl. In the case of everywhere positive potential 
U(r), this relation does not fulfilled definitely since, for 
such a potential, the quantities u n defined in i|25|) should 
be connected by a reciprocal inequality following from 
the Cauchy-Schwarz-Bunyakovsky inequality for the in- 
tegrals (|25p. Thus, using the degenerate kernel (|26|) in- 
stead of the true kernel U(k — k'), one must assume that 
the conditions (|29|l and l|35|) are fulfilled. 



VI. GENERAL SOLUTION OF THE 
SELF-CONSISTENCY EQUATION 

Substituting the kernel (|26[) into the equation (|14|l and 
grouping the terms independent of the momentum and 
also linear and squared terms one can see that the de- 
pendence of the order parameter on the momentum of 
the relative motion of the pair has the form similar to 
the Eq.® 



A(fc) 



( X k) - bk 2 



(36) 



where a and b are unknown parameters and x is an un ~ 
known constant vector being subjects to be determined. 



Thus, one can see that the solution of the equation 1|14|) 
reproduces the momentum dependence of the degener- 
ate kernel. Therefore, to solve this equation, one has 
to determine the parameters a, b and X- To find these 
parameters one has to substitute the kernel l|26fl and an 
explicit form of the solution (|36[) into the equation JT3J. 
After that, one obtains a system of integral equations 
determining the parameters a, b and %. One of these 
equations, following from the comparison of the coeffi- 
cients at the first power of k in the left and right sides of 
the equation resulting from the Ea. p4|l . has the form 



f{k)A(k)kd 2 k. 



(37) 



The equation l|37[) determines the vector x implicitly 
because the right-hand side of this equation depends on 
A which itself depends on X- It is n °t difficult to see 
that there is at most one solution of the equation 137|l . 
To prove this proposition let us introduce an auxiliary 
function H (A) = A • /(77(A)) which is a monotonically 
increasing function of its argument. Writing down the 
equation (|37l) in the component-wise form one can see 
that its right-hand side is a decreasing function of a com- 
ponent \i °f the vector x (i = 1, 2) whereas the left-hand 
side increases with \%- Therefore the equality l|37|l turns 
out to be possible at a single value of the vector param- 
eter X- 

The obvious solution of l|37|l is x = 0- Indeed, due to 
central symmetry of the domain S, the right-hand side 
of i|3 7|) vanishes at x = as a result of an integration of 
an odd function over a symmetric domain. It should be 
noted that the order parameter A(fc) and the excitation 
energy £(— k) = £(fc) are even functions of k at x = 0. 
There are no other solutions of the equation l|37() . Thus, 
in the case of a degenerate kernel Q26H , the equation 114|) 
may presuppose a simple solution of the form 



A(fc) =a-bk 2 , 



(38) 



determining the SC gap as a parabolic function of the 
momentum of the relative motion of the pair. 

Taking into account that x = one can write down 
two equations determining the parameters a and b: 



J a -\ji + l)a-(j l -^J 2 )b = u. 



1 



-J Q a 



1 



Ji - 1 & = 0, 



where 



f(k)k 2n d 2 k 



(39) 



(40) 



where n — 0, 1, 2. Thus, we have a quasi-linear system of 
the equations (|39|l which contains three positive integrals 
(|40|l depending on two unknown parameters a and b. 

First of all, let us multiply the inequality k 2 < 1/2 by 
the function f(rj(k)) and then by the another function, 
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k 2 f(r](k)). Integrating the obtained inequalities over the 
domain 2 one can deduce two new inequalities, J 1 < J /2 
and J 2 < J%/2>, which allow to conclude that both coeffi- 
cients in the first equation (|39l) are positive. Therefore, 
the solutions of the equation system (|39|l have one and 
the same sign (positive, for example). As a result, the 
SC order parameter as a function of the momentum of 
the relative motion changes its sign at a certain k — fc 
and can be written as 



A(fc) = b{kl - k 2 ) 



(41) 



where k^ = a/b. The parameter b determines an energy 
scale of the SC gap and k is a radius of the circle on 
which the SC gap vanishes. In the case of repulsive inter- 
action considered here, it is obvious that the SC gap has 
to vanish inside the domain S because at U(\k — k'\) > 
there exists no constant-sign non-trivial solution of the 
equation i|14|) . 

The conditions that a > and b > 0, as it follows 
from the second equation i|39|) . lead to a restrictive in- 
equality J 1 > 2 which is a consequence of more strong 
inequality following from the condition of the non-trivial 
consistency of the system (I39[) . Calculating the determi- 
nant of the system (|39|l one can write down this condition 
in the form 



{J 1 /2-lf = J (J 2 /4-l). 



(42) 



It is seen that it must be J 2 > 4. Therefore, taking into 
account that J l > 2J 2 , one can obtain J l > 8. 

Instead of the pair parameters a and b characterizing 
the SC gap it is convenient to consider another pair, b 
and k . To determine these two parameters one can use 
any pair of the equations (|39f) and (|42|l . As one of them 
it is convenient to use the second equation of the system 
(|39|l writing it in the form 



Jl 



k 2 J 



(43) 



Then, expressing the left-hand side of the equation l|42() 
with the help of the second of the equations (|39|) . the 
Eq.l|42p can be rewritten as 



Jo 



(44) 



The system of equations <|43[) and 144|) is fully equivalent 
to initial system (|39|l and can be used to determine the 
parameters b and k of the SC gap. In this connection, it 
should be noted that the integrals J , J 1 and J 2 , because 
of their dependence on A, depend themselves on b and 



VII. ORDER PARAMETER 

Let us rewrite the equations 143|) and l|44|) in the ex- 
plicit form 



1 



(k 2n - kl n )d 2 k 



47T J s ^2( fe ) +6 2( fc 2_j fc 2)2 



= 2 r ' 



(45) 



where n = 1,2. Is is impossible to solve the equation 
system (|45f) analytically in the whole of the domain of 
definition of unknown quantities (b > 0, < k\ < 1/2) 
because the integrals in l)45|) depend on the form of the 
domain of integration. However, it is easy to determine 
these integrals in the two limiting cases of large (b — » oo) 
and small (b — > +0) values of the gap energy scale b. In- 
deed, at large b one can neglect the augent in the radicand 
of each of the integrals (|45|l . After that, these integrals 
can be calculated easily and turn out to be proportional 
to 1/6 therefore one can obtain the solution of l|45|l in a 
trivial way. 

Considering the integrals in l|45|l in the most interesting 
case of small b one has to take into account the fact that 
pair excitation energy JQ) vanishes on the line separating 
occupied, and unoccupied, parts of the do- 

main S as far as this line (PFC) represents a piece of the 
FC. This results in a divergence of the integrals (|40ll at 
small b and, if one takes into account that the main con- 
tribution into the integrals is from a narrow strip along 
the line £(fc) = 0, it is this divergence that results in the 
existence of the solution of the system <|45ll . Let us de- 
note a width of the strip as 2AA; and, taking account of 
the equation of the line £(fc) = in polar coordinates, 
k = k(ip), consider an integral of a general form 



1 

J = — 



F(k) k dk dip 



(46) 



where F{k) is a continuous function without any sin- 
gularity. Let us consider a behavior of this integral at 
6^+0 using a strip of the width 2Afc as an integra- 
tion domain. First of all, one can perform an integration 
over a polar angle in the limits from ip 1 up to p 2 corre- 
sponding to the endpoints of the line £(fc, p) — 0. Then, 
introducing a new variable of integration £, one can per- 
form an integration over k in the limits from k(p) — Ak 
up to k(ip) + Ak. On account of the smallness of Afc, one 
can assume that k and £ are connected with each other 
in a linear way: £ ~ (d£(fc(</?)) /dk)(k — k{p)). In addi- 
tion, the argument k in the integrand can be replaced 
by its constant (at a given ip) value corresponding to the 
PFC: k = k((p). Thus, after integration over £, one can 
obtain a singular (logarithmic in b) contribution into the 
integral (|4*6*|l . 

Besides the singular contribution there is also a regular 
one which depends on one-particle dispersion and on the 
form and size of the domain S. This contribution depends 
on kg as on a parameter. Thus, the integral (|46|l can be 
represented as 



J = Aln(l/6) + C, 



(47) 



where the coefficient in front of the logarithm has the 
form 

1 r* 3 kF{k) 



-dip, 



(48) 



where an integration over the angle variable should be 
performed along the PFC (see Fig.l). The regular part 
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of the integral denoted as C can be taken equal to its 
value at 6 = 0. 

Thus, at small 6, the equation system 1(45(1 can be writ- 
ten in the form 



(A n -k^A )Hl/b) = P n 



(49) 



where the coefficients A n are independent of k 2 , and can 
be determined by the expression l|48(l in which one has 
to put F(k) = k 2n whereas the right-hand sides of the 
equations (|4*9")l depend on kg and can be written as 



P. 



h,2n (~i 



(50) 



Here C n are the regular parts of the integrals l(4t)(l. 

As it follows from the Eas. (|49() . the equation sys- 
tem 1(45(1 has a non-trivial solution with necessity. To 
prove such a statement let us consider each of the equa- 
tions H45f) as an equation of a line on the plane b, k^: 
b n = b n (k 2 ). Then the point of the intersection of these 
two lines corresponds to the solution of the equation sys- 
tem 1(45(1 . Putting fcp — in these equations one can 
obtain the left-hand sides of each of them as functions of 
b (Fig. 3). These functions monotonically decrease from 
+oo at b — > +0 to zero at b — > +oo. Due to the inequal- 
ity 1(25} . the second of these functions, corresponding to 
n = 2, with necessity is less than the first one (n = 1) 
at any b. Therefore, at k 2 = 0, the solution ^(0) of 
the first equation turns out to be more than the solution 
b 2 (0) of the second equation because of the fact that the 
right-hand side of the second equation is always more 
as compared to the right-hand side of the first one. In 
the another limiting case when b —> +0, one can use an 
asymptotic representation of the equations l|45|l in the 
form 1(49|> . As far as the right-hand sides of the equa- 
tions H49[) are finite and In (1/6) — > +oo at 6 — » +0 the 
values of fcg corresponding to 6 — > can be written as 
= A 1 /A at n = 1 and k 02 = \J A 2 j A a at n = 2. 
Taking into account the condition l(29|l one can see that 
&o n < 1/2. There is Cauchy-Schwarz-Bunyakovsky in- 
equality, A\ < A A 2 , for the quantities A , A x and A 2 
defining by such integrals as (|48|l therefore < k$ 2 and, 
with necessity, there exists a point of the intersection of 
the lines 6 ra = 6 n (/c§). 

To obtain the non-trivial solution in the case of small 
6 one has to perform term-by-term division of the equa- 
tions H49fl . This leads to a closed equation determining 
the parameter k$. As a result, a characteristic energy 
scale of the SC gap can be found from any equation of 
the system ((49(1 and, for example, can be written in the 
form 



exp 



Pi 



A 1 - A k Q 



(51) 



This relation explicitly determines the dependence of SC 
order parameter on the momentum of the relative motion 
of the pair. Due to the fact that 6 << 1, it is obvious 
that the exponent in l|51|) is negative. 




FIG. 3: Schematical illustration of the graphic solution of 
the equation system (1491 . 



One can see that, in the case of small 6, the line k — k Q 
intersects the PFC. Indeed, let us suppose that 6 and 
k Q are the solutions of the equation system (|45|) . At 
6 — > +0 the integrals in 1)45(1 are determined by their 
singular parts following from the integration over a small 
vicinity of the PFC. An assumption that the circle k = k 
does not intersect the PFC means that both the PFC and 
its small vicinity are situated either inside or outside this 
circle. Therefore in the whole of this vicinity one may 
expect either the inequality k < k or the opposite in- 
equality k > k . In such a case, the integrand in l|45|) 
must be of constant sign and the integral corresponding 
to n = 2 has to be less than the integral with n = 1 be- 
cause of the inequality k 2 + k 2 . < 1 following from 12911 . 
Thus, if one consider that the circle k = k and the PFC 
does not intersect each other the equations in the sys- 
tem (|45fl turn out to be incompatible with each other. 
Consequently, at 6 — > +0, the parameter k turns out to 
be such that the intersection of the PFC and the circle 
k = k arises with necessity. In such a case, the quantity 
k 2 — fcg has different signs on different pieces of the PFC 
interior and exterior the circle k = k . Such a conclusion 
proved under the condition that 6 — > +0 valid also at 
rather small but finite values of the parameter 6. 

As far as the circle k = k intersects the PFC, one can 
assume that there exist the minimal, k , and maximal, 



values of fc corresponding to the endpoints of the 
PFC. Taking into account the explicit expressions of the 
integrals l|48|l at F(k) = k 2n , one can obtain an estima- 
tion 



\A 1 — fcoA)l — 



Aip k\ 



M "'m 



(52) 



2tt (v F ) 

where ~ ip 2 — ip 1 is an angular size of the PFC and 
1_ J_ p dip 
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is an average (over the PFC) value of the dimensionless 
Fermi velocity. One can see from Ea. (|52ll that the abso- 
lute value of the SC order parameter decreases exponen- 
tially with a decrease both of the angular size and the 
anisotropy (characterizing by the difference k 



M 



J Of 

the PFC. In particular, in the case of repulsion-induced 
pairing with zero total pair momentum (when the PFC 
coincides with the FC), the non-trivial solution similar 
to Ea. l|51[) is absent. 

Thus, using conventional energy and momentum di- 
mensions, one can write down the SC order parameter, 
being the non-trivial solution of the self-consistency equa- 
tion l|14|). in the form 



A(fc) = A (fc)exp 



(54) 



where 

w = (A 1 - A Q kl)/P x (55) 

can be treated as an effective interaction strength and 

A Q (k) = U r 2 (k 2 ~k 2 ) (56) 

is a pre-exponential factor dependent on the momentum 
of the relative motion of the pair. 



VIII. STEP- WISE APPROXIMATION OF THE 
ORDER PARAMETER 

The continuous solution Eq. 1)54)1 of the self-consistency 
equation can be approximately presented as a step-wise 
function of the momentum of the relative motion^S^ 
Here, we briefly discuss the simplest form of such an av- 
eraged solution assuming that the true solution, Eq. (|54"jl , 
is replaced by an approximate one which has the form 
A(k) = A p , p = 1,2, where 



A(k)d 2 k. 



(57) 



The parameters A 1 and A 2 defined by the Ea. (|57(l may 
be considered as average values of the order parameter 
inside the domain S at k < k and k > k respectively. 

Performing the averaging of the left- and right-hand 
sides of the self-consistency equation 1)140 one can obtain 
the equation system 2 ^ 

2A 1 = -aHE^A^ - (1 - a)EU 12 A 2 f 2 , 

2A 2 = -aEU^AJ, - (1 - a)EU 22 A 2 f 2 . (58) 



where a = 5^/3, 1 — a = r E. 2 /E, and 

1 f d 2 k 



fp f V (^23 



(2n) 2 E p J 3p /£2 (fc) + A 2 



(59) 



is the value of the functional, 



A(k)d 2 k 



(60) 



corresponding to the average value A p of the order pa- 
rameter inside the subdomain 3 p . The average value of 
the interaction matrix element U(\k — k'\) corresponding 
to a momentum transfer k — k — k' is defined as 



Uppi Upip ^ 



U{n)d 2 K 



(61) 



where S* p , is a domain of definition of the momentum 
transfer k when an initial momentum, fc, belongs to the 
subdomain S p whereas a final one, k' , belongs to S / . 

If one represents the functions l|59l) in the explicit form, 
the equations 1)58)1 lead to a closed system of two tran- 
scendental equations. The integrals 1)59)1 over the sub- 
domains S p are similar to the integrals 1)40(1. Therefore, 
there exists a logarithmic singularity of such an integral 
when the domain of integration includes the PFC. The 
analysis of the non-trivial consistency of the Eqs.JSSJ) 
shows that the non-trivial solution exists under the con- 
dition that 



^12^21 " U U U 22 > 0. 



(62) 



This condition coincides with the well-known Suhl- 
Matthias- Walker condition — and is definitely fulfilled in 
the case of repulsive potentials which have an least one 
negative eigenvalue. 

The approximate approach developed here to deter- 
mine the order parameter arising at singlet repulsion- 
induced pairing and corresponding to extended s -wave 
symmetry can be used in the case of d -wave symmetry 
as well. In this connection, it should be noted that a 
simple replacement^ of the true kernel U(k — k') of the 
self-consistency equation by a constant 

U(k - k') -> {A(k)U(k - k')A(k'))/{A 2 (k)) (63) 

treated as an average over the FC hardly ever leads to a 
non-trivial solution because of the fact that such an av- 
eraging procedure, similar to performed in Eas. l|16|l and 
(I23|l , means the summation over both k and k' resulting 
in a non-negative value of the constant l|63|l in the case of 
repulsive interaction. Thus, strictly speaking, consider- 
ing the scattering due to a repulsive interaction U(k-k'), 
one has to take into account not only a difference be- 
tween the areas of the momentum space where the order 
parameter has different signs but also the dependence of 
the interaction matrix element on the momentum trans- 
fer k — k' . 



IX. CONCLUSION 

In this paper, we have shown that there must be at 
least one negative eigenvalue of the interaction matrix 
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U(k — k') to ensure the possibility of a repulsion-induced 
non-trivial solution of the self-consistency equation. A 
rise of the pair Fermi contour results in the fact that such 
a necessary condition of repulsion-induced superconduct- 
ing pairing becomes the sufficient one at arbitrarily weak 
interaction strength. 

It should be emphasized that negative eigenvalues have 
to arise with necessity in real Fermi systems, in particu- 
lar, due to Friedel oscillations of screened Coulomb po- 
tential. Quasi-two-dimensional electron system typical of 
cuprate superconductors may be considered as an inter- 
mediate case between comparatively weak Friedel oscil- 
lations in three-dimensional system and one-dimensional 
charge density wave behavior. The pair Fermi contour 
arises when, at certain large total pair momentum, there 
exists mirror nesting of the Fermi contour. 

Repulsion-induced singlet pairing results in an ax- 
tended s-wave symmetry of the superconducting order 
parameter. This parameter is defined inside the do- 



main of definition of pair relative motion momentum and 
changes its sign on a line intersecting the pair Fermi 
contour. An anisotropy is a necessary feature of one- 
electron dispersion which leads to the repulsion-induced 
non-trivial solution of the self-consistency equation: the 
corresponding superconducting order is absent in the 
framework of any isotropic model. 
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